Atom gratings produced by large angle atom beam splitters. 
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An asymptotic theory of atom scattering by large amplitude periodic potentials is developed in 
the Raman-Nath approximation. The atom grating profile arising after scattering is evaluated in the 
Fresnel zone for triangular, sinusoidal, magneto-optical, and bichromatic field potentials. It is shown 
that, owing to the scattering in these potentials, two groups of momentum states are produced rather 
than two distinct momentum components. The corresponding spatial density profile is calculated 
and found to differ significantly from a pure sinusoid. 
PACS number(s): 03.75, 03.75.D, 39.20, 8I.16.T 



I. INTRODUCTION 

Interaction of atoms with a spatially varying optical 
field can lead to a spatial modulation of the atomic den- 
sity. For example, when a coUimated atomic beam propa- 
gating along the x axis passes through a resonant, stand- 
ing wave field directed along the z axis, it acquires trans- 
verse momentum components 2nhk and {2n + l)hk for 
the atomic ground and excited states, respectively, where 
k is the field propagation constant and n is an integer In- 
terference between these momentum states results in 
a periodic modulation of the atom density with period 
d = A/2, where A = 27r/fc. This splitting of the atomic 
wave function (resonant Kapitza-Dirac effect |^ ) and the 
density modulation have been observed in a beam of Na 
1^ and Ar jj] atoms, respectively. Atom gratings can be 
deposited on substrates and have potential applications 
as diffraction gratings for soft x-rays. 

The grating density profile p (z) produced by this tech- 
nique contains an infinite number of Fourier harmonics 
having periods A/2n and amplitudes p„. For the grating 
profile to be approximately sinusoidal, then at most a few 
Pn contribute in the Fourier sum. On the other hand, the 
"sharp" grating that can be produced by atom focusing 
using strong standing wave fields consists of a large 
number of Fourier components. 



X/w > 1, 



(1) 



where w <C A is the spot size of the focussed atoms. 
Sharp gratings have important lithographic applications, 
but are not well suited for use as elements to diffract x- 
rays. For x-ray scattering, one finds that each density 
component p„ leads to scattering at angle 



±nXx/d, 



(2) 



where A^, A is the x-ray wavelength. If one moni- 
tors scattering at a specific angle, most of the scattered 
energy does not contribute to the signal. As a conse- 
quence, it is most efficient to scatter x-rays from a sinu- 
soidal atom grating, since only one Fourier component 



is nonvanishing [rt = 1 in Eq. (|^)] . To increase the an- 
gle and resolution of scattering one needs a high order 
sinusoidal grating having period 



A/2n, 



(3) 



where n > 1. Several techniques have been proposed to 
produce such high order gratings. One can distinguished 
two groups of methods, harmonic suppression (HS) and 
large angle beam splitters (LABS) . In HS, one arrives 
at a periodicity (||) as a result of suppressing all Fourier 
components in the atom wave function that are not a 
multiple of n. Included in the HS group are methods that 
employ photon echoes , multicolor fields Q , high order 
Bragg scattering and Fourier optics techniques [ p^ , 
methods which are summarized below. 

An appropriate echo geometry for HS involves atoms 
that pass through two standing wave fields, spatially sep- 
arated by a distance L (or atoms that are subjected to 
two temporally separated standing wave field pulses). At 
distances L' = [n' > n, where n' and n are integers 
having no common factors), the atomic density contains 
only those Fourier components having n that is a mul- 
tiple of n. The other components are suppressed if the 
atom beam has a transverse velocity width (56' > A/i 0. 
A ^—period grating {n = 2) has been observed using this 
method on *^Rb atoms cooled in a MOT and subjected 
to the 2 time-separated standing wave field pulses . In 
the multicolor field technique ||], atoms interact during 
time T with a traveling wave field having frequency il and 
two counterpropagating traveling wave fields having fre- 
quencies Q, + 5i and f7 + (52- If I^^tI ^ 1 and the detunings 
are chosen such as niSi + n252 = 0, then a grating having 
period A/2f;, is produced, with n = ni ri2. In contrast 
to the echo-technique, the multicolor field method allows 
one to obtain near-unity contrast for arbitrary h. For 
Bragg scattering, a proper choice of the angle between 
the atomic beam and the field propagation direction or, 
alternatively, a proper choice of atom-field detunings , 
leads to situations where low spatial harmonics are sup- 
pressed. A -g-— period grating has been obtained using 
this technique In the Fourier optics technique one 
focuses atoms transmitted through a standing wave field 
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using a long-focus lens. Different harmonics of the atomic 
wave function are focused at the different positions. Us- 
ing spatial filtering, one can suppress all harmonics ex- 
cept two, whose interference leads to the grating of the 
desirable period. In this way one achieves 100% contrast 
pure sinusoidal gratings, but the amplitude of the grat- 
ings decreases with increasing grating order. 

An ideal large angle beam splitter (LABS) splits an 
initial state having transverse momentum p = into two 
states having momenta 



p — ±nhk, 



(4) 



where n ^ 1. The interference between these states leads 
to a spatial density 



p (z) = 1 + cos {2nkz) 



(5) 



having period X/2n. The prototype LABS consists of a 
V-potential 



U{z,t) = fit) hkn\z\, 



(6) 



off of which atoms scatter. The function /(i), normal- 
ized such that /^^^2/(0'^^ — describes the potential 
experienced by an atom in its rest-frame [t — x/u, where 
u is the projection of the atom velocity along the beam 
propagation direction [x axis) and r is the pulse duration 
in the atomic rest frame]. In the Raman-Nath approxi- 
mation 



(7) 



where 



is a recoil frequency and m is the atomic mass, one can 
show that, at a time t following the atom field interac- 
tion and for n 3> (wkt) , an atom grating is produced 
having density 



hk 

p{z) ~ 2 [1 -f cos (2nfcz)] for |z| < n — t; 

^ 1 for \z\ > n — t. 

m 



(8) 



In the past several years, a number of methods have 
been proposed for using counterpropagating optical fields 
to create potentials that simulate a V-potential (|^). 
Atom scattering off these fields involves multiphoton pro- 
cesses that result in a transfer of momentum between the 
counterpropagating field modes. The effective potential 
U {z,t) corresponding to these processes is a function of 
bi-linear products of the counterpropagating field am- 
plitudes. As a consequence, U {z,t) has period A/2 as 
a function of z. One might expect that a so-called tri- 
angular potential, a periodic potential coinciding with a 
V-potential within a period. 



U{z + X/2,t) = U{z,t), 



(9) 



should produce gratings similar to the sinusoidal grat- 
ing (H) produced by the V-potential. Based on this idea 
several schemes for LABS have been proposed. 

The simplest example of such a beam splitter is the 
strong standing wave field first considered in ||]. After 
adiabatic elimination of the internal degrees of freedom, 
one finds a sinusoidal potential having period -j. Asymp- 
totically, for a large potential amplitude, the main con- 
tribution to the scattering arises from the area around 
points zo where the potential slope is maximum; in this 
region the potential is approximately linear. The momen- 
tum an atom acquires after scattering near such points is 
p ~ Ft, where F ^ ~ i^U/dz)^^^^ is the classical force 
acting on atom. 

An alternative method which produces a potential that 
more closely approximates (||) is the magneto- optical field 
geometry [ p^ , in which resonant cross-polarized counter- 
propagating fields drive a ground-to-excited state transi- 
tion, J = J — 1, where J is the atomic angular mo- 
mentum. A magnetic field directed along the z axis that 
splits the excited state Zeeman sublevels is also applied. 
A regime has been found where the magneto-optical po- 
tential is extremely close to the triangular potential (||). 
The nonlinearity near z — zq is small compared with that 
associated with the standing wave potential. Magneto- 
0£tical forces have been used to defiect an atomic beam 
|l3| and to split an atomic beam into two momenta com- 
ponents localized near values p — ±21?ifc [n — 21 in Eq. 
®] [0- Atom scattering from a magneto-optical po- 
tential is considered in ||l5|-[l7|], while the force acting on 
atoms is calculated in fig] . Instead of using a magnetic 
field one can detune standing waves symmetrically with 
respect to the atomic transition frequency . 

Another way of creating a potential similar to the tri- 
angular one is to use a bichromatic field in which 
atoms scatter from two standing wave fields, shifted from 
one another with respect to their frequencies, and spa- 
tial and temporal phases. Forces resulting from this field 
have been observed in [^,|2j . A remarkable property of 
these forces (large magnitude in comparison with radi- 
ation forces, insensitivity to atomic velocity) allows for 
fast atomic beam deceleration |^ and for large angle 
deflections ||2^ . The triangular potential associated with 
a bichromatic field was predicted in |25|, and scattering 
off this field has been observed |26| and studied for two- 
level and three level atoms p8|-^. Moreover, it 
has been shown that atoms can be scattered with 
unit probability with a given change of momentum if a 
"counterintuitive pulse sequence" is used. 



2 



Despite the similarity between the triangular poten- 
tial (||) and V-potential , there is an important differ- 
ence. Owing to the intrinsic A/2 periodicity, atoms scat- 
ter in the potential (||) not only into states p = ±nhk but 
also into states having p = ± (n ± 2) Tik, ± (n ± 4) ?ifc .... 
Thus, at best, the potential (^) can produce two groups 
of momenta centered at p = ±nhk. As such the resulting 
pattern deviates from that of a true V-potential. The 
additional momentum states produce deviations from an 
ideal sinusoidal atom grating by producing additional 
Fourier components. 



Piy (z) = fli/ cos [2 (n + v) kz] 



(10) 



where ^ 1 for v ~ ±An/2 and An ^ n is the range of 
components created by the potential. Therefore, LABS 
produces atom gratings that are not pure, high order 
gratings; their overall periodicity is still A/2. Gratings 
arising in HS also contain extra terms (p^, but one is 
guaranteed that <C ao if is not a multiple of n. As 
a result, atom gratings produced by HS closely approxi- 
mate high order ^-period gratings, having an envelope 
that oscillates slightly with period A/2 around a constant 
value. 

Scattering into a range of momentum components oc- 
curs even for the triangular potential. For potentials 
which approximate the triangular potential (standing 
wave, magneto-optical scheme, bichromatic field), addi- 
tional features arise owing to the nonlinearities in these 
potentials. It turns out that the width An for these po- 
tentials increases with increasing n. If one tries to de- 
crease the grating period by increasing n, there is a corre- 
sponding deterioration of the degree to which the grating 



profile approximates a pure sinusoidal grating. 

In this article we provide a detailed calculation of 
the spatial gratings produced by LABS, using analyti- 
cal asymptotic and numerical calculations. We consider 
only the Fresnel regime when all scattered atomic beams 
overlap with one another. An asymptotic expression for 
the atomic wave function in momentum space is derived 
in the next Section and the atomic density profile is cal- 
culated in Sec. 3. Results are summarized and discussed 
in Sec. 4. In the Appendix, we consider atom gratings 
produced by the triangular potential (S). 



II. ATOM WAVE FUNCTION 

In this section, we calculate the atomic response to a 
standing wave field pulse, a magneto-optical field pulse 
and a bichromatic field pulse, all in the Raman-Nath ap- 
proximation. For an atomic beam, the pulsed interaction 
occurs in the atom rest frame. 



A. Standing wave. 

When a two-level atom interacts with a standing wave 
field pulse having electric field. 



E {t) = yE^ [t] cos {kz) e 



+ c.c. 



(11) 



having (real) amplitude EQ{t), propagation constant k, 
duration r, and frequency w, the atomic wave function 



(in the interaction representation) evolves as 





X (t) COS {kz) exp (—iAt) 



X (t) cos {kz) exp (iAt) \ 

j ' 



(12) 



where ip and V'e are ground {g) and excited (e) internal states wave functions, A = w — is the atom-field detuning, 
X {t)=iJ,Eo (t) /h is a Rabi frequency, is the ground to excited state transition frequency, and ^ is a (real) dipole 
moment operator matrix element associated with the e ^ g transition. In Eq. ( p^ ) we neglect a term, {jp' /2mfi) '9, 
[p is the atomic center-of-mass momentum operator), which is valid in the Raman-Nath approximation (|^). For an 
adiabatic pulse 



by *_ = 
the wave function is '4'j_ = 



x/x<max{|A|,|x|}, (13) 
^ . If the wave function before the pulse arrives is given 
then, immediately following the pulse, the atom returns adiabatically to its ground state and 



it is convenient to use a semiclassical dressed state approach 
V- [z) ■ 


[z) 





where 



V'+ {z) 
viz) 

e{z) 




cos {2kz)) 



1/2 



(14a) 
(14b) 

(14c) 



T] {z) is a transmission function and 9 (z) is the spatially inhomogeneous pulse area. To simplify the calculations, 
it is assumed that the pulses considered in this paper have leading and trailing edges of duration Tr, and that the 
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pulse Rabi frequency is constant and equal to x between the leading and trailing edges. Adiabaticity is guaranteed if 
|A| Tr ^ 1. Moreover, if t^/t <C 1, the contribution to the phase 9 {z) from the tails is smaller by a factor of order 
Tr/r from that of the central region and modifies only slightly the points of stationary phase to be calculated below. 
As such we approximate 9 (z) as 

9{z) = 9w{z), (15a) 
w{z) = {1 + e [1 + cos (2fcz)]}^/^ (15b) 

where 

9 = At/2 

and 

i = 2xVA'- 

If initially an atom is in the zero momentum state, then the Fourier components 

i>n= -TTT exp [-2inkz - i9{z)] (16) 

Jo -^/^ 

determine the amplitudes for scattering into momentum states p — 2nhk, which we refer to as scattering into state n. 
In the asymptotic limit 

^max ^min » 1, (17) 

the main contribution to the integral (|l^ ) arises from points of stationary phase, z — zq. These contributions are 
maximal if the second derivative of the phase also vanishes at z = zq. This means that the dominant contribution to 
scattering into state n occurs for n « no , where no is determined by the equations 

no = -9' (zo) /2fc, (18a) 
9" (zq) = 0. (18b) 

Expanding the phase in ( p^ to third order in (z — zo) , one finds that ipn is given by the asymptotic expression, 

V'n ~ 2fc [2/ \9"' (zo)|]^^^ exp [-2inkzQ - i9 {zq)] Ai {An2A: [2/9"' (zo)]^^^} , (19) 



where 

An = n — no 

and Ai (x) is an Airy function. For the field parameters 
given in (Iq), one has 



An, 



no = ±no; 

no-0{[l + e-(l + 2O'/'] /2}'^', 
9"' = ±8k^no, 

zo = ±Acos-i{[(i + 2e)-i-e] /a 

9{zo) ^9{l + 20'^\ 
such that 

\4>n\' - At' [An(2/no)'/'' 



(20a) 
(20b) 
(20c) 
(20d) 
(20e) 

(21) 



The distribution ( ^l| ) is plotted in Fig. |^ and has width 
(full- width at half maximum) 



1.29n, 



1/3 



(22) 



The result (^ij) differs from that obtained in reference . 

Using an asymptotic expression for the Airy function 
for large negative arguments one finds that for An < 
and I An I ^ 1 the distribution oscillates as a function of 
An. Averaging over these oscillations one finds 



|^„|' ~ vr-i (2noAn)-^/^ 



(23) 



which coincides with the momentum distribution associ- 
ated with classical particles scattered by a standing wave 
potential ||||. 

In the case of a far detuned field (^ ^ 1) , when 



>{z) =0|l + |[l + cos(2fcz)] j>, 
no = 9^/2, 



(24a) 
(24b) 
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one finds from Eq. (|6|) that the exact expression for the 
Fourier harmonics is 



i^n = Jn m/2) exp [-imr/2 - i9{l + C/2)] 



(25) 



where J„ (x) is a Bessel function. The asymptotic re- 
sult (|l9|) follows directly from Eq. (^5|) using a well- 
known asymptotic expression for large order Bessel func- 
tions One can consider Eq. (19) as a generalization 
of the asymptotic expression for Bessel functions in the 
limit that n ^ 1 for a function having arbitrarily large 
periodical phase. 



B. Magneto-optical potential 

This regime of scattering arises |T3| , when atoms in- 
teract with an optical field 



-iujt 



+ C.C. 



(26) 



and a static magnetic field along z. If the optical field 
drives aJ = 0— !-J=l ground-to-excited state transi- 

tion, then the wave function '3/'= tp-e , ("0' — ipe^^*', 

i/j-e, i^e are amplitudes for the ground state, mj = — 1 
excited state, and nij — 1 excited state, respectively) 
evolves as 



where the interaction Hamiltonian is 

/A G- G+ 
KU ^hlGl -LUL 
\G1 LUL 



(27) 



(28) 



zbo-iL are Zeeman shifts of the states ±e, 

G± = T2"^/^Wfle±"/4 cos (fcz T 7r/4) , 

ujfi = 2 X 3~^/^i?o {t) is a Rabi frequency, and is a 
reduced matrix element of the dipole moment operator. 
To develop an asymptotic solution of Eq. ( p7| ) , one must 
first solve the characteristic equation 



\U -Ml = 0. 



(29) 



Solutions can be found in Refs. ||15| , |16yi8| . Let us denote 
by A (z, t) the root of Eq. ( p9| ) which equals A when the 
field vanishes. If initially an atom is in its ground state 
and if there is no level crossing during the interaction, 
then after the pulse the atom returns to its ground state. 
The ground state wave function ip following the pulse is 
given by equations analogous to Eqs. (|lj) with a pulse 
area defined by 



(z) = / dt[X{z,t) - A] 



(30) 



A potential that closely resembles the triangular poten- 
tial (^) emerges if one chooses a rectangular pulse having 





uiR — 2ujl and A = 



(31) 



(although A = 0, adiabaticity is still maintained if 
i^LTr S> 1). For these values, we find from Eqs. ( |2^ ) 
and (ISOl) that 



(a) 




(b) 




(c) 



60 



'40 



■20 



20 



40 



60 



p/hk 



FIG. 1. Probabilities l-^nl for atom scattering by LABS 
into momentum states p = nhk located at n ~ ±no for 
no = 207r/3. corresponding to the modulation amplitude cho- 
sen in ||l5| for numerical calculations; (a) standing wave field, 
(b) magneto-optical field, (c) bichromatic field. Solid lines are 
the numerical values and the gray solid lines are asymptotic 
expressions obtained from Eqs. ^ |42| ) . 
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6{z) = -0w{z), (32a) 
w (z) = sin I i sin~^ [sin (2fcz)] | , (32b) 

where 9 = 2ujlt. The function w [z) has discontinuous 
derivatives. In the regions of continuity this function is 
given by 



w (z) — sin ( 2^-2^ j : for < z < — , 



sm 



sm 



for - < z < -A, 

,3^ A 
' 8^<"< 2- 



(33a) 
(33b) 
(33c) 



The slope of w (z) has extrema at zq = and zq = A/4, 
associated with sphtting of the incident state in two mo- 
mentum groups located at no = ±noj with 

no = 61/3. (34) 

From Eq. (|l|) it follows that 

~ (±1)" (18/no)^^^ Ai [±An (IS/no)^^^! . (35) 



Exact and asymptotic values for \'4'n\ as a function of n 
are shown in Fig. [|. One sees that asymptotically, for a 
given angle of splitting (i. e. a given value of no), the net 
result of replacing a standing wave field by a magneto- 
optical field is to increase the maximum of the scattering 
probability by a factor S'*/^ and decrease of the distribu- 
tion width An„io by a factor 3^/^. In spite of this, the 
width of distribution is still an increasing function of the 
angle of scattering, 



0.62n, 



1/3 



(36) 



C. Bichromatic field 

Finally, we consider the interaction of a two-level atom 
with a quantized field composed of two standing wave 
fields having frequencies cjo ± A, the same initial num- 
ber of photons V \, and shifted with respect to one 
another by a quarter of a wavelength. The Hamiltonian 
can be written as [ETI 



H — Ho + Hi, 

Ho = hujo (|e) (e| -I- c^c -I- dU) , 

Hi ^ hA (c^c - d^d) - fig {cos (fcz + tt/A) [c^ \g) (e| + |e) (5I c] -i- cos (fcz - 7r/4) [d^ \g) (e| 



(5|d]} 



(37a) 
(37b) 
(37c) 



where c and d are annihilation operators for modes 
± A, (7 is a coupling constant, and Hi is the inter- 
action Hamiltonian. If initially an atom is in its ground 
state, then in the interaction representation the system 
wave function is given by 



/ 



ipg {1^ - + l;z, t) 
■0e [v - 1, v; z,t) 
ij;{z,t) 

Ipe {v, V -l\ Z,t) 

ipg {v + - l;z,t) 



\ 



(38) 



where (m, n; z, t) is the amplitude for the atom to be 
in state \i) and the field to be in a state having m and 
n photons in the modes ujq± A, respectively. Before the 
atom field interaction, ■0 {z, t) = -tpg {v, v; z, t). This wave 
function evolves as 



(39a) 



2A 

A 



n 



+ 







-A 



\ 



0+ 
-2A 



(39b) 



where n± = Qq cos (fcz ± 7r/4) , and CIq — g^pu is a Rabi 
frequency. Let us denote by A (z, t) the eigenvalue of the 
matrix L associated with a dressed state (eigenvector of 
the matrix L) which coincides with the initial state of the 
system before and after the interaction. As in the cases 
above, if adiabaticity is maintained, the initial state wave 
function undergoes a phase shift as a result of the inter- 
action. The final state wave function is still given by Eq. 
(|l4b| ), with Q (z) = dt\ (z, t) . 

A remarkable situation had been found in reference 
H. With a Rabi frequency f^o = V3A/2, (z) closely 
approximates a triangular potential. Explicitly, one finds 



where 



Q (z) = ~Qw (z) , (40) 
-At and the function w (z) is shown in Fig. 
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0.125 



0.375 



0.5 



0.25 

FIG. 2. One period of the normalized pulse area 
w (z) = 6 (z) /6 for scattering by a bichromatic field. 

The slope of w (z) has extrema at 2; = and z — A/4, 



where ' 



±0.624fc, 



±0.211fc'^. Using these values 



one finds that atoms scatter in two states near ho 
with 

no w 0.3126*, 



The corresponding Fourier components are 
V-n (±1)" 2.870 rig ^/^Ai (=F2.870An/no/^ 

In this case the distribution |V-'n|^ has width 
Arifcc = 0.57riy^. 



±no, 



(41) 



(42) 



(43) 



This distribution is plotted in Fig. |^. The distribution, 
obtained here by a numerical integration of Eq. (|l^) with 
pulse area (|40), differs from the distribution calculated in 
the article p7| by a numerical solution of the Shrodinger 
equation. We suppose that differences arise owing to the 
fact that a rectangular pulse profile was used in [ p7| ; such 
a pulse does not satisfy our adiabaticity requirements. 



III. ATOM GRATING PROFILE 

The LABS split an initial p = state into two groups 
of states having momenta close to p = ±2nohk. It is 
of interest to compare the matter gratings generated by 
these LABS in the Fresnel scattering zone with those of 
an ideal beam splitter. An ideal beam splitter produces 
purely sinusoidal gratings having period 



At a time t after scattering from a potential having 
period A/2, each momentum state's amplitude evolves 
as V'n exp (— iw2fc?T-^i) , where Lj2k — 2hk'^ /m is a recoil 
frequency. The atom density p{z,t) is given by 



p (z, t) ~ ^ V'nV'n' exp \2ikz {n — n') — iu>2kt {r 



(45) 



Asymptotic calculation of the atomic density profile is 
simplest for a triangular potential (|^) (see Appendix). 
One finds in that case that the grating profiles are not 
those of an ideal beam splitter. If one wishes to carry 
out analogous solutions for the potentials considered in 
Sec. II, purely numerical calculations are needed. Such 
calculations are not given in this paper. 

We are left with the problem of evaluating Eq. ( |45| ) 
in some asymptotic limit. Immediately we run into two 
problems: (1) Owing to the nonlinearities in the poten- 
tials, one cannot simply use the asymptotic expression 
for the Fourier components (^^ since, when these ex- 
pressions are substituted into Eq. (45), one finds that 
the constraint. 



(46) 



is violated. Equation ( p^ follows from Eq. (^5|) and 
the fact that 2 J^^^ p {z,t) dz/ \ = 1. If one restricts the 
summation to values n Ki ho, the series X^IV'nl^ calcu- 
lated using Eq. ( ^5| ) diverges as J2 ~ "0)^^^^ ■ This 
example demonstrates that an alternative technique is 
needed to calculate the atomic density in the asymptotic 
limit, which allows one to include contributions from all 
momentum states. (2) Potentials having smooth rather 
than triangular extrema lead to the atom focusing 
at a time t ~ l/9uj2k |36j, i. e. to a sharp l-period grat- 
ing having little in common with a sinusoidal grating. 
There is no simple way to predict the times for which 
this density pattern is transformed into one that might 
resemble that of the triangular potential. Rather than 
looking for such times numerically, we consider an alter- 
native approach for obtaining the grating profile, valid 
for large times. 

To this point we assumed implicitly that we were deal- 
ing with a monovelocity atomic beam. The time t is 
related to the observation plane at a distance x from the 
interaction zone, hy t = x/u, where u is the atom ve- 
locity along the a;-axis. In principle, for a LABS one 
needs a monovelocity beam, since the pulse area 9 and 
(therefore) the angle of scattering (no) depends on u like 
1/m; however, one can neglect this dependence for a small 
width of the velocity distribution, Au ^ u. Although we 
require Am <C u, it is assumed that Au is nonvanishing 
such that, for 



di = A/4no. 



(44) 



(47) 
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the grating is an oscillating function of u \ ^T\ . 

On averaging over u one finds that the phase factor 
exp \—iu!2kt (n'^ — n''^)] appearing in p5| ) averages to 
zero unless |3q| 



or antisymmetric pulse area [9 (z) = ±6 (—z)] , one can 
simplify Eqs. (49c, 49d) as 



*W ^l[F{z) + F{X/4-z)], 



(48) 



F(z) 



Gratings of this type have been observed recently in the 
optical Talbot effect 

Setting n' — ±n in Eq. (^) one finds that the overall 
grating period becomes A/4, instead of A/2. Inside a pe- 
riod of the grating profile (0 < z < A/4) , one arrives ||3^ ] 
at the following expression for the density; 



dz 



exp [idx [z, z)] , 



A/4 

x{z, z) = w (z + z) — w {z — z) , 



(50a) 

(50b) 
(50c) 



p{z,t) 
B 

$(z) 
F{z) 



B + $(z), 

^'^ dz_ 
A/2 



h(^)l'- 



A/2 



[F {z) + F{z + A/4)] 



dz 



z\<Min{z,^-z} •^/4 



r/ (z + z) rj* {z — z) . 



(49a) 
(49b) 

(49c) 
(49d) 



where the functions w (z) characterizing the various po- 
tentials have been given in Sec. II. 



A. Standing wave field 

To simplify the calculations, we consider the asymp- 
totic grating profile for a far-detuned standing wave 
field = 2x^/A-^ ^ l) whose pulse area is given by Eq. 

( P4a| ) . For this pulse w{z ) = {l + |[l + cos(2/cz)]} and 



X (z, z) = S_ sin {2kz) sin (2fcz) 



(51) 



where _B is a constant background of order unity and 77 (z) 
is the transmission function defined in Eq. (14a). We will 
refer to $ (z) as the grating profile. Although $ (z) can be 
negative, the atom density p (z, t) is always positive. For 
a phase transmission function (14b) having symmetric 



For l arge no — 0^/2 ^ 1, one can evaluate the integral 
( ^Ob| ) using a stationary phase method. The points of 
stationary phase are at z = ±A/8. These points belong 
to the region of integration if z > A/8. Thus, neglecting 
terms of the order of 1/9, one finds 



F (z) ~ e (z - A/8) 2^/2 ]^^Q^ sijj (2fcz)]"^/^ cos [6*^ sin (2/fcz) - 7r/4] 



(52) 



where (z — A/8) is a Heaviside step-function. This expression becomes invalid near z = A/4; however, for A/4 — z ^ 
A/no the integral ( ^Ob| ) can be approximated as 



F (z) - 2 Jo [9i (tt - 2fcz)] 



Substituting Eqs. ( |52| , |53| ) in (50a) one obtains 

$ (z) '-^ Jq (4nofcz) , for z < A/no; 

^ [ttuq sin (2fcz)] ' cos [2no sin (2fcz) — 7r/4] , for z 7^ and z 7^ A/4; 
- Jo [2no (tt - 2A:z)] , for A/4 - z < A/no. 



(53) 

(54a) 

(54b) 
(54c) 



The profile <I> (z) is shown in Fig. ||. For large no one sees from ( |54b| ) that $ (z) oscillates rapidly with a spatially 
inhomogeneous period equal to 



d (z) = di/ cos (2fcz) . 



(55) 



Although there are rapid oscillations with reduced periodicity, the grating does not approach the purity associated 
with an ideal beam splitter. 



B. Magneto-optical potential. 

For the magneto-optical scheme the pulse function of z, has discontinuous derivatives at multiples of 

z = A/8 [see Eqs. ( ^ ^3|)]. Since both arguments z ± z of the w-functions in Eq. (50c) vary in the interval [0, 2z] for 
Izl < z, one has to consider 4 cases: 



a 



) z e [0, A/16] , ^) z e [A/16, A/8] , 7) z g [A/8, 3A/16] , S) z £ [3A/16, A/4] . (56) 
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0.125 0.25 
Z/A. 

FIG. 3. Grating profile produced by a far-detuned standing wave fie ld. Solid gray lines are the exact profile calculated 
using Eqs. ( |50| , dashed and dot-dashed lines are asymptotic results ( |54a , 54c) and (54b), respectively. 



For the simplest case a) one uses expression (33a) for both functions w {z ± z) to find x {z, z) — 2 c os ^k z sin |fc 



Stationary phase points occur for Zg = ±|A ^ [— z, z], and one can use an integration by parts in Eq. (50b) to arrive 
at the asymptotic expression 



F{z) - 3 sin 



' sm I —kz 
3 



nd cos^ ( —kz 



(57) 



We will see belovif that only for small z ~ l/k9 and for \z — A/4| ~ l/kO will terms of order 9^^ be important. In all 
other cases terms of order unity (arising for values of z where x {z, z) ^ 1/6*), and terms of order 9^^/^ (arising from 
integration near stationary phase points) dominate. If such terms are absent, we say that corresponding contribution 
vanishes. 

In case (3) the function w {z + z) is equal to sin (z -|- z)] and sin [i (tt — 2k {z + z))] in the intervals z G = 



[- 



■zA 



z\ and z G i?+ = — ^, z] respectively, while the function w (z ~ z) is equal to sin [|/c {z 



and 



sin (tt — 2k{z — z))] in the intervals z £ A- = [z — ^,z'\ and z e = [—2, 2 — f] . Since C A± one has to 
map out the region of integration [— z, z] into 3 intervals: 



(58) 



Let us denote by Fi (2) the contribution to F (z) from interval i (i = a, b, c) For z G a, x (z, z) 
2 sin [i (4fcz — tt)] cos (tt -|- 4fcz)] , z^ = —A/8 ^ a, but the function x becomes small near z = A/8, i. e. 





X' 




A A 




'A 


a = 




, 6 = 


z , 2 

8' 8 


, c = 





r/6 



A 



— 2 cos ( 



Fa (z) ^ — I d(p exp 
0, otherwise. 

The functions Fi, (2) = 0, (2) = F* (2) , and one finds for case f3) 



, for - - 2 < 1/ke, 
8 



(59a) 
(59b) 
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"4 /A 
F{z) ^jy0ki^--z 

~ 0, otherwise, 



for - - z < 1/9, 



where 



6 r^^ 

J {z) — — I c(0cos|zcos( 



(60a) 
(60b) 

(61) 



Function (61) is plotted in Fig. ^. 



IN 




FIG. 4. J(z) ( |6l| ) as a function of z near the points 2 = (n + ^) 



In case 7) one has to distinguish intervals a — [— z, -I — z] , 6 = [-1 — z, z — -1] , c = [z — -I, z] . When z <E a, 
X (z, z) = 2 sin (4fcz — tt)] cos [| (4fcz + tt)] , Zg — —A/8 G a and the leading contribution to Fa {z) is 



Fa (^) 3 <^ ttB sin 



(4fcz - tt) 



-1/2 



exp < 2i0 sin 



(4A:z - tt) 



while Fb (z) = 0, Fc (z) = F* [z) . Adding these contributions one finds for case 7) 
F{z) ~2X(z), 



K{z) =i{TT 



{Akz - tt) 



-1/2 



cos I 26 



{Akz - tt) 



(62) 

(63a) 
(63b) 



In the region near A/8 [0 < z — A/8 < 1/^] , where Eq. ( p^ ) is invalid, one returns back to expression (60a). 
Another feature arises in case S), when one has to map out the integration region z e [— z, z] into 5 intervals: 





3 ■ 






a = 


-z,z--A 




z 



2- ttA, ■ 



A 




'A A' 




A 


z 

8 


, c = 


z, z 

8 8 


, d = 


z 

8 



-A-z 



-A — z, z 



(64) 



The stationary phase point Zg belongs to the region of integration for the contribution Fb (z) , where x ^) = 
2 sin [g (Akz — tt)] cos (4A:z + tt)] , Zg = —A/8 and F^ (z) is given by the right-hand-side of Eq. ( |6^ ) . To achieve 
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this result one needs the region of integration to be larger than X/O^^'^ ^ but the interval h in Eq. ( |6^ ) has zero length 
at z = A/4, rendering the method of stationary phase invalid. In this region the function Fb (z) is simply equal to the 
integrand a,t z = Zg multiplied by the length of interval b : 



Fb (z) w ^kSz exp i9 (1-2 ■ 3"^/2fc(5z) , for z = A/4 - Sz, Sz < X/O. 



(65) 



Calculating all other contributions using integration by parts, one arrives at the following final result for region near 
z=A/4 



F{z) = — i 2kSzecos 
no ' 



1 - 2 • 3-^/^kSz 



sin -OkS 



Substituting F (z) into Eq. (jSOal) one finds [0 < z < A/8] 
2 

$(z) - J 



$(z) - £^ |2fcz6'cos 



4 

-ekSz 
3 



(1 - 2 • 3-^/^ki 
X 



1/2 > 



(1 - 2 • 3-i/2fcz) } , for z < X/e; 



, for z = Sz, Sz < X/0; 

8 

$ (z) ^ K (z) otherwise. 



(66) 

(67a) 

(67b) 
(67c) 



where the functions J (a;) and K (z) are given by Eqs. (|6l| 63b ) respectively. The profile for A/8 < z < A/4 can be 
obtained using the symmetry property, $ (|- — z) = $ (z) . The profile ( |6^ ) is shown in Fig. ^. 
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As in the case of the standing wave beam sphtter, everywhere except near the points z = 0, A/8, A/4, the profile is 
a grating having a spatially inhomogeneous period 



d (z) — di j cos 

The profile is not that of an ideal beam splitter. 

C. Bichromatic beam splitter 

Lacking an analytical expression for the potential pro- 
duced by a bichromatic field we are not able to obtain an 
analytical asymptotic expression for the grating profile. 
On the other hand, it is possible to evaluate Eqs. (50a) 
numerically. In Fig. ^the profile is shown for rip — 207r/3 
and is similar to that for the magneto-optical potential. 



e 




p — ±2nohk, 



(69) 



- {Akz - tt) 

D 



to a pure, sinusoidal high-order grating 

p (z) = 1 + cos (Anokz) 

having period 



(68) 



0.125 0.25 

z/A. 

FIG. 6. Density profile produced by bichromatic field for 
no = 207r/3. 



IV. DISCUSSION 

Asymptotic expressions in the limit of large-angle scat- 
tering have been obtained for the atomic wave function 
and grating profile after scattering from a periodic trian- 
gular potential (see Appendix) and from potentials which 
approximate this triangular potential. The potentials are 
chosen |l^j2^ with the goal of producing two and only 
two final state momenta 



di = X/Auq. 



(70) 



(71) 



with (riQ S> 1) . In turn, interference between these mo- 
menta components in the Fresnel scattering region lead 



In practice, it is impossible to produce only two fi- 
nal state momenta with periodic potentials. At best, 
one finds scattering into two groups of states centered 
around p = ±2no?ifc. Owing to the range of momenta in 
each group, none of the potentials can generate the ideal 
grating profile (|70|). For the potentials considered in this 
work, the width of the momentum groups scales asymp- 
totically as Hq^^ [see Eqs. (^, 43)]. If one's goal is to 
creates a higher order sinusoidal grating, then one needs 
to increase the angle of splitting, i. e. the parameter uq. 
However, by increasing uq one broadens the width of the 
momentum groups, and diminishes degree to which the 
grating profile approximates a pure sinusoidal grating. 

In the limit of strong fields and large scattering angles, 
we developed asymptotic techniques for evaluating the 
atomic wave function and density. One might think that 
a quasiclassical or WKB method could be used in this 
limit, but the linear nature of the potential invalidates 
such an approach in the region of maximum scattering. 
Instead, the method of stationary phase was used to ob- 
tain asymptotic expressions for the wave function. The 
results depend solely on the first and third derivatives of 
the potential at the point of maximum slope. For all but 
the triangular potential, however, these asymptotic ex- 
pressions cannot be used in the expressions for the atomic 
density since they do not produce a convergent result. 

Only for the triangular potential does the asymptotic 
expression for the Fourier amplitudes converge rapidly 
enough as a function of Fourier index n to allow one to 
calculate the atomic density using these asymptotic for- 
mulas. The periodic triangular potential is analyzed in 
the Appendix. The scattered signal is a periodic function 
of time with period to — T/2nQ, where T = 27r/cj2A; is the 
Talbot time (recall that time t is related to the distance 
X following the atom-field interaction zone hy t = x/u, 
where u is the longitudinal velocity). In a given scatter- 
ing plane following the atom-field interaction, the atom 
density consists of three distinct spatial regions, one of 
the form ( [70| ) but with twice as large an amplitude, one 
with p (z) = 1, and one with p (z) = 0. Examples of these 
profiles are shown in Fig. |^ in the Appendix. 
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For the standing wave, magneto-optical, and bichro- 
matic field potentials, the atomic density was calculated 
for scattering distances x 3> Tu'^/Au Effl, where Am <C u 
is the range of longitudinal velocities . The density is 
a periodic function of z having period A/4. 

For the standing wave potential one can distinguish 
regions about z = nA/4 (n is an arbitrary integer) of ex- 
tent X/np where the density has a sharp maximum [see 
Eq. (|54a| )] of order unity. With the exception of these 
regions the density is an oscillating function of z with am- 
plitude of order Uq and spatially inhomogeneous pe- 
riod ([55|). The period becomes infinite at z — {n + ^) j, 
which implies that, near these points, the profile is flat. 
All these features are seen in Fig. |^, where the grating 
profile is plotted for uq = 2 and no — 207r/3. The last 
case corresponds to the value chosen in for numeri- 
cal calculations. It is surprising that the analysis of the 
scattering in terms of two distinct spatial intervals seems 
to be valid even for no = 2, when the condition no 3> 1 
is not satisfied. 

For the magneto-optical case one has to distinguish 
three regions 



I I 
II. 



(n+i) 

-nA < 



< 



A/no, 



A/no, 

III. other values of z 



(72) 



accounts for most of the range when no ^ 1, the grat- 
ing profile is similar to that found for the standing wave 
potential, i. e. it is a grating having amplitude of order 

— 1/2 

Uq and a spatially inhomogeneous amplitude period 
given by Eq. (^8|). In contrast to the standing wave case, 
this period cannot be larger than 2 x 3~^^^di. The grating 
profile produced by the magneto-optical field is shown in 
the Fig. |. For a moderate value no = 2, the dependence 
( |67c) ) associated with region III is a poor approximation 
to the profile; however, one can see from Fig. ^ that 
it is not necessary to use the expression for region III. 
Regions I and II overlap with region III and with each 



other, and Eqs. (37b, 37a) provide a good approximation 
to the profile for all z. 

The grating profile for the bichromatic potential is 
shown in Fig. ^ for no — 207r/3. 

In summary, we have asymptotic expressions for the 
momentum distribution and calculated the grating pro- 
file for several classes of optical beam splitters. While 
the momentum distribution resulting from these beam 
splitters can approach that of an ideal beam splitter, the 
associated grating profiles do not closely approximate a 
pure sinusoidal pattern. 
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APPENDIX A: APPENDIX - TRIANGULAR POTENTIAL. 

It is instructive to consider scattering by the ideal periodic triangular potential (|^). This potential differs from 
those considered in the text in that the second derivative in each half-period of the potential vanishes identically. The 
pulse area associated with this potential is given by 

where no — n/2. For the Fourier components (|l^), one finds 

, _ r2nosin[(n - no)7r/2] /vr (n^ - ng) , for n 7^ -no 

'^""\(-l)"V2, forn = -no ' ^ ^ 

Consider now the case of large integer no. For n close to inp one finds 

"ipno+An ~ sin [An7r/2] /ttAh , (A3a) 

V'-«o+A„ ~ {-1}"" [S (An) - V„o+An] . (A3b) 



where S (An) is a Kronecker symbol. In a contrast to the nonlinear potentials considered above, the amplitudes (A3) 
decay sufficiently fast to calculate both the atomic wave function and atom density. In the expression for the wave 
function, 
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V' {z, t) = ij+ {z, t) + V- {z, t) , 

V'± {z, *) = X] V'±no+An exp 2ikz (±no + An) - iw2fe (±no + An)^ t 



An 



(A4a) 
(A4b) 



one can neglect terms quadratic in An in the exponents if |t — mT\ <C T, where m is an integer and T = 27r/a;2fe 
is the Talbot time. For these values of t, one can obtain analytic asymptotic expressions for the wave function and 
density. The wave function is a periodic function of z having period A/2 and a periodic function of t having a reduced 
Talbot period 



to = 7r/a;2/c?^o• 



(A5) 



Within one period 



ip+ {z, t) = exp {2inokz - iuj2knot) Q 
ip- {z, t) = exp {-2inokz — iu)2knot) © 



„ fik 
2nn — t 





„ hk 




( 


z + 2no — t 


-1) 




m 





(A6a) 
(A6b) 



where Q (x) is a Heaviside step-function. 

One sees that after scattering the wave function splits in two components having average momentum ±2nQhk. 
Moreover, there is a spatial modulation of the wave function corresponding to j sections of the incident plane wave 
scattered with momenta ±2no^. These wave function components, having unit amplitude and spatial width j, are 
centered at z = z±, where 



„ hk 
z+ = 2no — t - 
m 



+ 2' 



„ hk 
Z- = — 2no — t + 
m 



(A7a) 
(A7b) 



(n± are arbitrary integers). Interference occurs when different components centered at z± overlap. Outside of the 
interference regions, the density can equal unity for points belonging to one of the 2;± -regions or for other points. 
Analyzing all these possibilities, one arrives at the following density profile 



p {z, t) =^ {z, t) + 2 (-1)"" cos (Anokz) (j) {z, t) , 



(p{z,t) = < 



Uz,t) = 



3A 

'1. 



if z e 

^ a z G 

if z e 

if z e 

in other cases 

0, if z G 

1, if e 

2, if z G 



ta 
2 



- Zt, I + zt] and t G [--_ 
Zt,-^-Zt] andtG [-\,0] 
^-zt,^+zt] andie [0,f] 
-^ + zt,f-zt] andte 



0, if z G 

1, if z G 

2, if z e 

0, if z G 

1, if z G 

2, if z G 

0, if z G 

1, if z G 

2, if z G 



A _ 



.. 3A 
Zl, — 

3 A 



Zt] U [| + zt,-|-zt] U [f+zt, 



Zt, 



g + Z(, 



zt\ 



-i-t+^t]u[-|-z*,|+zt]u[| 



+ zt,-k - Zt\ 



Zt 



- - - Zt - - 

-|, -I - z'] U [-1 +zt,l- Zt] U [| + Zt, A] 

A A 



Zt, ■ 



Zt 



'X + t ~ Zt] 

1' tx'"^ U [| - Zt, -I + Zt] U[f- Zt, 
"S + zt,^ - Zt] 



for t G 



for t G 



for t G 



to 

' 2 ' 



t_o 

' 4 



for t G 



to_ io 
4 ' 2 



(A8a) 



(A8b) 



(ASc) 



(ASd) 



(A8e) 



(A8f) 



where 



Zt = 



Xt^ 
2to 
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This profile is shown in Fig. The function (j) (z, t) represents an envelope function for higher order gratings having 
period di — X/Auq. 




- 0.25 - 0.125 0.125 0.25 - 0.25 - 0.125 0.125 0.25 



z/X z/X 

FIG. 7. Scattering by the triangular potential. The first seven graphs give the density profile for a single period at different 
times. The last graph if the density averaged over to. The time to is the Talbot time divided by no and no = 11 in this figure. 



Expression (A8b) gives the envelope at a given time. One can consider also the grating profile as a function of t for 
given z. This would correspond to atoms arriving at a given z with different longitudinal velocities, allowing one to 
make a comparison with the results of Sec. III. For fixed z as a function of t, one finds 



(t>{z,t) 



if t e 
if t e 
if t e 
if t G 



-|to-2fto,ito + 2fto] andzG[-|,-|] 
-ito + 2fto,-iio-2fio] andze[-|,0] 
Ito - 2f io, ito + 2f to] and z € [O, |] 
-ito + 2fto,ito-2fto] andze[|,f] 



(A9) 



in other cases 



One sees that over one period of Talbot oscillations the high harmonic amplitude is constant for some times and 
vanishes for others. The harmonic amplitude averaged over a Talbot oscillation period is simply given by (0 (z, t))^ = 
tz/to, where tz is the duration of the time interval when the harmonic amplitude is constant; the background term 
on average is equal to unity [(^ (z, t))^ = 1] . Calculating tz from A9, one finds 



{p (z, t))^ =1 + 2 (-1)"° cos (Anokz) 

1 - i |z| , for \z\ > I 



for 



z < 



(AlO) 



This result could be obtained using Eqs. (ttOL pOf) of Sec. Ill with the pulse area (Al) 
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